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Abstract 

We  represent  the  modelling  of  a  thermoelectric  cooler,  which  is  designed  by  using  micromachining  and  thin  film  technology.  The 
cooler  fabrication  is  compatible  with  standard  semiconductor  technology.  Therefore,  it  can  be  integrated  in  microelectronic  circuits.  The 
most  important  parameters  of  the  device  like  cooling  power,  maximum  temperature  difference  and  optimum  current  density  are 
calculated.  By  using  thermoelectric  thin  films  with  high  efficiency  and  very  thin  SiC/Si3N4-membranes,  a  cooling  power  of  a  few 
milliwatts  or  maximum  temperature  difference  of  30-50  K  can  be  achieved.  ©  1999  Elsevier  Science  S.A.  All  rights  reserved. 
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1.  Introduction 

Thermoelectric  coolers  (Peltier  devices)  are  widely  em¬ 
ployed  in  microelectronics  to  stabilize  the  temperature  of 
laser  diodes,  to  cool  infrared  detectors  and  charge-coupled 
devices  (CCD),  and  to  reduce  unwanted  noise  of  integrated 
circuits. 

A  conventional  thermoelectric  cooler  usually  consists  of 
a  number  of  n-  and  p-type  bulk  semiconductor  thermoele¬ 
ments  connected  electrically  in  series  by  metallizations  and 
sandwiched  between  two  electrically  insulating  but  ther¬ 
mally  conducting  ceramic  plates  [1].  The  dimensions  of 
commercially  available  coolers  varies  from  about  4X4X3 
to  around  50  X  50  X  5  mm3.  Although  in  principle,  the 
dimensions  can  be  reduced  further,  the  preparation  of 
conventional  thermoelectric  coolers  is  a  bulk  technology 
and  is  incompatible  with  microelectronics  fabrication  pro¬ 
cesses. 

Thin  film  thermoelectric  devices  with  very  small  di¬ 
mensions  have  been  fabricated  using  microelectronics 
technology  as  generators  [2]  and  sensors  for  radiation  [3,4], 
electrical  AC  and  DC  power  [5-8],  flow  [9]  and  low 
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vacuum  [10,11].  However,  to  date,  such  thermoelectric 
devices  have  not  been  used  successfully  as  thermoelectric 
coolers  because  the  substrate  is  thick  compared  to  the  thin 
film  thermoelements.  Consequently,  the  cooling  achieved 
by  the  device,  when  operated  in  the  Peltier  mode,  is 
drastically  reduced  by  the  substrate’s  thermal  bypass  [12— 
14], 

Recently,  we  have  developed  very  thin  SiC/Si3N4- 
membranes  (200  nm  thick)  with  compensated  mechanical 
stress  on  bulk  silicon  wafers.  These  membranes  are  pre¬ 
pared  by  PECVD  of  the  SiC/Si3N4  films  and  by  mi¬ 
cromechanical  anisotropic  etching  of  the  silicon  wafer 
from  the  backside.  The  small  thermal  bypass  of  these 
membranes  facilitates  the  fabrication  of  a  thin  film  thermo¬ 
electric  cooler  with  sufficient  cooling  power  and  tempera¬ 
ture  difference.  We  report  on  the  modelling  of  such  micro¬ 
electromechanical  coolers. 


2.  Thermal  model  of  the  thermoelectric  cooler 

Fig.  1  shows  the  top  view  and  cross-sectional  view  of 
the  thermal  model.  The  cooler  consists  of  a  number  n  of 
thermocouples  connected  in  series  (thermopile)  and  ar¬ 
ranged  with  their  cold  junctions  (temperature  Tc)  around 
the  cooled  area  F ,  whereas  their  hot  junctions  are  in  good 
thermal  contact  with  the  heat  sink.  The  thermopile  and  the 
cooled  area  are  supported  by  a  very  thin  membrane  with 


reserved. 


0924-4247/99/$  -  see  front  matter  ©  1999  Elsevier  Science  S.A.  All  rights 
PH:  S0924-4247(99)00002-3 


96 


F.  Volklein  et  al. /Sensors  and  Actuators  75  (1999)  95-101 


(a) 


low  thermal  conductivity  in  order  to  achieve  good  thermal 
insulation  of  the  cold  junctions.  The  quadratic  membrane 
results  from  the  anisotropic  etching  of  the  silicon  wafer. 
The  bulk  silicon  rim  acts  as  a  heat  sink.  Heat  is  pumped 
laterally  from  the  central  region  to  the  silicon  substrate  rim 
and  then  dissipated  vertically  through  it  to  an  external  heat 
sink.  The  cooled  area  F  is  modelled  as  a  quadratic  area 
with  the  length  a  of  the  sides.  Each  thermocouple  is  a 
stack  of  two  thermoelectric  films  (A  and  B,  respectively) 
with  an  insulating  layer  between  them  (see  Fig.  lb).  The 
width  of  a  thermocouple  is  b,  its  length  (the  distance 
between  heat  sink  and  cooled  area  F)  is  l  and  the  spacing 
between  thermocouples  is  s. 

The  modelling  of  the  cooler  can  be  carried  out,  in 
principle,  by  three-dimensional  simulation  tools.  However, 
for  the  thin  film  device  of  Fig.  1,  the  analysis  can  be 
performed  in  terms  of  a  one-dimensional  model.  Usually, 
the  thickness  of  the  thermoelectric  films,  of  the  membrane 
and  the  insulating  film  is  in  the  order  of  1  |xtn,  whereas  the 
thermocouple  length  is  in  the  order  of  1  mm.  Therefore, 
temperature  distributions  in  the  thickness  dimension  can  be 
ignored.  Furthermore,  because  of  the  symmetric  design  of 
the  cooler,  only  one  separated  thermocouple  has  to  be 


analysed.  Such  a  thermocouple  represents  a  cantilever 
beam  with  a  cooled  area  F/n  at  the  tip  of  the  beam.  In 
such  cantilever  beams,  the  heat  expansion  can  be  described 
by  only  one  coordinate. 

Fig.  2  shows  the  cross-sectional  view  of  the  thermal 
model  for  such  a  single  thin  film  thermocouple.  We  as¬ 
sume  that  the  temperature  Tc  is  uniform  over  the  whole 
area  F/n.  This  can  be  achieved  by  a  metallization  of  F, 
which  also  reduces  the  absorption  of  thermal  radiation  and 
therefore  minimizes  the  impact  of  heat  from  the  ambient. 
Furthermore,  heat  is  absorbed  by  the  area  F/n  from  the 
ambient  due  to  conduction/convection  of  the  surrounding 
gas.  N  represents  an  additional  thermal  load,  for  example 
by  a  microelectronic  device  which  has  to  be  cooled.  The 
cold  junction  at  the  position  x  =  l  generates  the  cooling 
power  Q  (Peltier  ‘heat’),  which  is  given  by  the  Seebeck- 
coefficients  of  the  legs  (aA  and  aB,  respectively),  the 
current  I  and  the  temperature  Tc 

Q  =  (aA-aB)ITc.  (1) 

In  the  volume  element  dV  of  the  cantilever  beam  heat 
is  absorbed  by  radiation  and  convection  from  the  ambient, 
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Joule’s  heat  is  generated  by  the  current  I  through  the 
thermocouple  and  heat  is  conducted  from  the  heat  sink 
with  ambient  temperature  7a  to  the  cold  junction  with 
temperature  Tc.  The  heat  balance  equation  for  the  volume 
element  dV  at  position  x  with  a  temperature  Tc  <  T(x)  < 


where  yF  =  8  eF  a  7a3  +  2/1  denotes  the  heat  transfer  coef¬ 
ficient  for  radiation  and  convection  and  sF  denotes  the 
emissivity  of  the  area  F.  The  solution  of  Eq.  (2)  is 


\tdtb-Jp  cosh (/y^p  ) 


—  (8scrT*  +  2/3)7>(T(  x)  —  Ta)  =  0.  (2) 

Here,  the  first  term  describes  the  heat  conduction 
through  the  cantilever  beam,  the  second  term  is  Joule’s 
heat  and  the  third  term  is  the  heat  absorbed  by  radiation 
and  convection  with  a  heat  transfer  coefficient  y  =  8  scr  T2 
+  2/3.  We  assume  that  ( T(x )  —  7j)/Ta  <s:  1  and  we  use 
the  approximation  scr(T(x)4  —  Ta4)  =  4saT2(T(x)  —  Ta) 
for  the  radiation  term,  a  is  the  Stefan-Boltzmann  con¬ 
stant,  s  is  the  emissivity  of  the  beam,  /3  is  the  heat 
transfer  coefficient  for  convection  of  the  surrounding  gas, 
pA  and  pB  are  the  electrical  resistivities,  Aa  and  AB  are 
the  thermal  conductivities  and  dA  and  dH  are  the  thick¬ 
nesses  of  the  thermoelectric  legs.  The  term  A tdt  describes 
the  total  thermal  conductance  of  the  cantilever  according 


where  Am,  dm  and  A;,  di  denotes  the  thermal  conductivity 
and  thickness  of  the  membrane  and  insulating  film,  respec¬ 
tively.  The  differential  Eq.  (2)  can  be  solved  by  using  the 
boundary  conditions 


(T(x)  -Ta)\x=0  =  0 
d(T(x)-7;)| 


=  ~{Q-N)-yF- 


I2R  ■  exp(  — /y^p  ) 

“  /^cosh(Vp)  Sinh(X^) 

+  7^(1-exp(-x^))  (6) 

where  R  denotes  the  resistance  of  the  thermocouple  and 
p  =  y/(A tdt).  Using  Eq.  (6)  for  the  calculation  of  T(x  = 
l)  -  Ta  =  7'a  -  7a,  we  find 


\  cosh(/y^p )  I 

=  -(e-w)+/fiW+4vFT'  (7) 

I  \,dtb-Jp  F\ 

Here,  the  factor  G  =  - — 1=—  +  yF—  of  (7’a  —  7j) 

I  tanh  [l\p  J  n  I 

represents  the  thermal  conductance  of  the  thin  film  device. 


The  factor  C  =  - 1= - 7 — i=z —  of  Joule’s  heat  repre  - 

/y  p  tanh(/y p  ) 

sents  the  exchange  of  Joule’s  heat  with  the  ambient  and 
with  the  cold  junction.  With  these  designations,  the  heat 
balance  equation 

(Tc  —  ra)  •  G  =  -  (  Q  -  N)  +  I2R  ■  C 

=  —  (  aA  —  aB  )  ITC  +  N  +  I2R  •  C 


(5) 


(8) 
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corresponds  to  the  equation  for  a  bulk  thermoelectric 
cooler.  However,  for  bulk  coolers,  the  thermal  conductance 
is  given  by  G  =  (Aa//Wa)  +  (AB//WB),  where  WA  and 
WB  are  the  cross-sections  of  the  thermoelectric  legs  and 
C  =  1/2. 

Using  Eq.  (8),  we  deduce  the  optimum  current  /opt, 
which  maximizes  the  temperature  difference  Tc  —  Td  and 
the  cooling  power: 


2RC 


(9) 


With  /  we  calculate  the  maximum  temperature  dif¬ 
ference 


(rc-ra)max-G  = 


( 


~«B  fTc2 

ARC 


+  N. 


(10) 


Without  any  heat  load  N,  the  thin  film  device  can 
achieve  a  maximum  temperature  difference 


(  Tc  -  Ta)max 


(«A-«B)2  rc2 

47?GC 


(11) 


The  maximum  heat  load  can  be  removed  from  the 
cooler  for  (Tc  —  Ta)max  =  0.  Consequently,  we  deduce  from 
Eq.  (10): 


ATmax 


(aA-  atB)2T* 
ARC 


(12) 


3.  Results  and  discussion 

We  discuss  the  maximum  temperature  difference  (Tc  — 
ra)max  and  the  maximum  heat  load  ZVmax  for  various 
operation  conditions. 

3.1.  Minimum  heat  exchange  by  convection  (operation  in 
high  vacuum) 


Then,  we  can  find  a  practicable  leg  length  Z,  for  which  the 
term  ( l-Jp  )  <K  l .  Since 


ybl 

\dtb/f? 


XtcLtb/l 


(13) 


the  condition  (ItJ p)  1  corresponds  to  a  small  thermal 
conductance  by  radiation  in  comparison  with  the  thermal 
conductance  by  heat  conduction.  For  e  =  1  (worst  case!), 
Ta  =  293  K,  dA  =  dB=  1  p.m,  ds  =  1  |xm  and  \tdt  = 
A AdA  +  A BdB  +  A sds  =  4.56  X  10“6  W/K  we  calcu¬ 
late  ( l-Jp )  =  /(1.58/mm).  Here,  A sds  =  A mdm  +  \Idi  de¬ 
notes  the  substrate’s  thermal  conductivity  and  thickness, 
including  the  thermal  conductivities  and  thicknesses  of 
membrane  and  insulating  film.  Typical  material  properties, 
which  we  have  used  for  calculations  are  represented  in 
Table  1.  For  /  <  0.3  mm  and,  consequently,  (l-Jp)  <  0.5 
the  approximations 


os  h(/yV) 


...  m2 


and  tanh  (l-\fp)  ~  (/\/p  ) 


can  be  applied.  With  these  approximations,  the  thermal 
conductance  is 


A  ,d.b  F 

G=  — —  +  yF- 


(14) 


the  factor  C  =  1  /2,  and  the  maximum  temperature  differ¬ 
ence  can  be  expressed  by 

(Tc  ~  /.Lax 


(aA  -  aB)2Tc2 


A  tdtb  F\ 

~r  +  r%) 


A  AdA  +  XBdB  +  A  sds  +  yF 


(15) 


For  operation  in  high  vacuum  (pressure  <  10  3  mbar) 
the  heat  flux  by  convection  can  be  neglected  ( /3  =  0). 


In  order  to  calculate  the  maximum  temperature  differ¬ 
ence,  we  have  to  study  the  thermoelectric  properties  of  thin 


Table  1 

Thermoelectric  properties  and  film  parameters  used  for  calculations 


Materials  or  geometrical  parameters 

at  Tc  =  250  K 

at  Ta  =  293  K 

Seebeck-coefficient  aA  =  —  aB 

=  a  (|xV/K) 

185 

230 

Electrical  resistivity  pA  =  pB  =  / 

o  (to-5  n  m) 

1.33 

1.70 

Thermal  conductivity  Aa  =  AB  = 

=  A  (W/m  K) 

1.28 

1.08 

Thermoelectric  efficiency  zA  =  t 

:B=Z  =  a2/(pAXlO-3/K) 

2.00 

2.88 

Emissivity  s,  sF  (worst  case!) 

1 

1 

Substrate’s  thermal  conductivity 

As  (W/m  K) 

2.0 

2.2 

Length  of  the  legs  l  (mm) 

0.3 

0.3 

Lateral  length  of  the  cooled  area 

a  (mm) 

1 

1 
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ds  =  1  pm 
ds  =  2  pm 
ds  =  3  pm 
ds  =  4  pm 
ds  =  5  pm 


0  2  4  6  8  10 

leg  thickness  d  (pm) 

Fig.  3.  Calculated  maximum  temperature  difference  as  a  function  of  leg  thickness  d  and  substrate  thickness  ds. 


film  materials,  p-type  (Bi0  5Sb1 5)  Te3-films  show  the  best 
thermoelectric  efficiency  at  room  temperature.  They  have 
been  prepared  by  flash  evaporation  [15]  or  MBE  [16]  and 
subsequent  annealing.  The  properties  of  these  optimized 
films  are  close  to  the  bulk  material  and  are  represented  in 
Table  1 .  The  preparation  of  n-type  films  with  a  high  figure 
of  merit  is  much  more  complicated.  Our  modelling  is 
based  on  the  (somewhat  optimistic)  assumption  of  a  high 
efficient  n-type  film  with  the  same  properties  like  the  best 
p-type  films.  By  further  improvement  of  the  preparation 
techniques,  certainly  it  will  be  possible  to  realize  such  high 
efficient  n-type  films.  Furthermore,  the  calculations  are 
performed  for  identical  leg  thickness  dA  =  dB  =  d.  Then, 
the  maximum  temperature  difference  can  be  expressed  as 

V — ^ — Try  (16) 

2r  +  2  Ad  +  y"2nbAdl 

The  leg  width  b  is  limited  by  the  circumference  V  =  A  a 
of  the  cooled  area.  For  a  very  small  spacing  s  the 
maximum  leg  width  is  given  by  bmax  =  U/n.  Using  bmax 
and  F/U  =  a/4,  we  learn  that  (Tc  —  Ta)max  is  independent 
on  b: 


(  Tc  -  Ta)m. 


(17) 


Fig.  3  shows  the  calculated  maximum  temperature  dif¬ 
ference  as  a  function  of  the  leg  thickness  d  and  substrate 
thickness  ds ,  calculated  according  to  Eq.  (17)  for  Ta  =  293 
K  with  the  film  data  at  Tc  =  250  K  (temperature  of  the 
cold  junction).  It  should  be  mentioned  that  in  Eq.  (17),  the 
contact  resistance  of  the  thermocouples  is  neglected. 
Therefore,  the  temperature  difference,  which  can  be  ob¬ 
tained  in  practice  will  be  smaller  than  the  calculated  one. 
The  curves  of  Fig.  3  represents  the  ultimate  limits  for  thin 
film  cooler  with  respect  to  the  recent  facilities  of  thermo¬ 
electric  films.  On  the  other  hand,  the  contact  resistance  of 


such  photolithographically  patterned  thin  film  devices  is 
only  a  small  part  of  the  total  resistance  R.  Therefore,  the 
contact  resistance  will  not  drastically  reduce  the  calculated 
curves. 

Using  Eq.  (12)  and  the  data  of  Table  1,  we  deduce  the 
maximum  heat  load  for  a  single  thermocouple: 


4  a2T2bd  a2T2bd 

2(Pa  +  Pb)1  Pl 


(18) 


and  the  total  heat  load,  removed  by 

4a2T2da 

N=n ■  Nma,  = - . 

pl 


thermocouples. 


(19) 


For  given  dimensions  of  the  cooled  area  (a  =  1  mm) 
and  of  the  leg  length  (/  =  0.3  mm),  N  can  be  increased  by 
increasing  leg  thickness  d.  This  is  due  to  an  increase  of  the 
optimum  current  /opt  with  increasing  leg  thickness.  The 
optimum  current  and  optimum  current  density,  respec¬ 
tively,  are  given  by 


7opt  ( aA  — aB)T  aT 
7opt  ~~bd~  2RCbd  ~  ~fA 


(20) 


where  we  apply  the  temperature  T=  Tc  for  the  case  of 
maximum  temperature  difference  and  T  =  Ta  for  the  case 
of  maximum  heat  load.  Using  the  data  of  Table  1  at 
Ta  =  293  K,  the  maximum  heat  load  and  the  corresponding 
current  density  of  the  thin  film  cooler  are  calculated  and 
represented  in  Table  2. 


Table  2 

Maximum  heat  load  N  and  optimum  current  density  jopt  of  the  thin  film 
cooler  as  a  function  of  leg  thickness  d 

Leg  thickness  d  (pm)  Heat  load  N  (mW)  y'opt  (103  A/cm2) 

1  L56  L32 

5  17.8  1.32 

10  35.6  1.32 
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Table  3 

(r„  —  Ta)max  and  maximum  heat  load  N  as  a  function  of  leg  thickness  d  for  a  leg  length  /  =  1  mm,  substrate  thickness  ds  =  1  p,m  and  with  y=yF=  100 
W/m2  K 


d(  p-m) 

Ufp) 

tanh  (l-fp) 

1 

<TC  ~  Ta)max  (K) 

./opt  <A°2  A/cm2) 

N  (mW) 

jopt  (102  A/cm2) 

cosh  (iffy 

1 

4.68 

1.000 

0.981 

-10.3 

9.20 

2.50 

9.28 

2 

3.75 

0.999 

0.953 

-14.4 

7.27 

4.01 

7.43 

3 

3.21 

0.997 

0.919 

-16.8 

6.15 

5.15 

6.36 

4 

2.86 

0.994 

0.886 

-18.4 

5.46 

6.11 

5.67 

5 

2.60 

0.989 

0.852 

-19.6 

4.94 

6.95 

5.15 

3.2.  Maximum  heat  flux  by  convection  (operation  at  nor-  calculated  with  the  parameters  of  Table  1  at  T  =  Tc.  Fur- 
mal  pressure  of  the  surrounding  gas )  thermore,  the  maximum  heat  load  (for  Tc—  Ta  =  0) 


By  using  the  microcooler  in  a  normal  atmosphere,  heat 
flux  by  convection  has  to  be  involved.  Usually  a  micro¬ 
cooler  chip  is  attached  on  a  header  (see  Fig.  1).  Applying 
4-in.  silicon  wafers,  the  distance  g  between  membrane  and 
header  is  about  500  p,m.  Then,  the  heat  transfer  coefficient 
1 6  for  convective  heat  flux  between  membrane  and  header 
can  be  expressed  by  /3=  A gas/g,  where  Agas  denotes  the 
thermal  conductivity  of  the  surrounding  gas  ((3=  52 
W/m2  K  for  air).  Furthermore,  free  convection  on  the 
other  side  of  the  cantilever  beam  and  thermal  radiation  has 
to  be  involved.  Therefore,  y=  yF~  100  W /m2  K  is  a 
reasonable  value  for  operation  at  normal  atmosphere.  Us¬ 
ing  this  value,  we  calculate  the  terms  (l^p),  tanh  (hfp ) 
1 

and  (1 - for  a  leg  length  /  =  1  mm  (Table 

cosh(/ypJ 

3). 

We  learn  from  Table  3  that  below  a  leg  thickness  of  4 
p.m  the  term  (l-/p)>  3  and  the  terms  tanh  (l-Jp)  and 
1 

(1 - — ,  ...•)  can  be  approximated  by  unity.  There- 

coshf/ypj 

fore,  we  deduce  from  Eq.  (7)  with  the  designations  of 
Section  3.1 

(Tc  -  TJ ( \tdtb y7  +  y/)  =  -(£?-A0+  l2R~(Jjpj 

(21) 

and  with  y=yF  for  the  maximum  temperature  difference 


(  Tc  -  7'a)max 


4 


1  + 


(22) 


Table  3  shows  also  the  maximum  temperature  differ¬ 
ence  (for  ds  =  1  p,m)  and  the  corresponding  optimum 
current  density 


.  4pt  (*a-<*b)T  <*T{p 

Jopt  bd  2RCbd  2  p 


(23) 


2a2Tfad\fp 
N=n ■  Amax  = - — 
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(24) 


and  the  corresponding  optimum  current  density,  calculated 
with  the  parameters  of  Table  1  at  T=Ta,  are  shown  in 
Table  3.  The  results  demonstrate  that  a  microcooler  using 
thermoelectric  films  with  z  =  2.00  X  10_3/K  can  perform 
a  maximum  temperature  difference  of  about  20  K  at 
normal  atmosphere.  Therefore,  e.g.,  dew-point  sensors 
based  on  such  microelectromechanical  coolers  can  be  real¬ 
ized  by  batch-process  fabrication. 
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